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Foreword 


G.  I.  Taylor  found  a  solution  of  the  gas  d;^maniic 
equations  which  describes  the  pressure  waves  produced  by  an 
explosion.   This  description  can  be  used  at  distances  from 
the  explosive  large  compared  to  the  explosive  dimensions, 
but  small  enough  that  the  shock  pressure  is  large  compared 
to  atmospheric  pressure.   In  fact  the  explosive  is  assumed 
to  be  a  point,  and  the  atmospheric  pressure  is  neglected 
compared  to  the  shock  pressure,  in  obtaining  the  solution. 
For  these  reasons,  the  solution  is  called  either  the  "point 
blast"  or  ''strong  shock"  solution.   The  same  solution  was 
also  obtained  by  J.  von  Neuinann,  and  a  solution  of  the  same 
type  was  found  for  explosions  in  water  by  H.  Primakoff . 
Similar  solutions  were  studied  by  G.  Guderley  in  Germany.' 

The  two  opposing  restrictions  on  the  range  of  validity 
of  the  Taylor  "point  blast"  solution  are  such  that,  for 
ordinary  explosives,  there  is  practically  no  range  in  which 
they  are  both  satisfied.   For  nuclear  explosives,  however, 
the  size  of  the  explosive  is  so  small  that  there  is  a  range 
in  which  this  solution  is  useful.   This  has  been  demonstrated 
by  comparing  this  solution  with  experimental  results. 

In  order  to  extend  the  solution  to  greater  ranges, 
where  the  shock  is  weaker,  it  is  necessary  to  take  account 
of  the  atmospheric  pressure  ahead  of  the  shock.   The  present 
report  by  Dr.  Morawetz  attempts  to  do  this  by  determining  a 
correction  to  the  "point  blast"  solution,  which  correction 
is  due  to  the  previously  neglected  atmospheric  pressure. 
A  system  of  linear  equations  is  obtained  for  the  correction 
and  their  solution  is  analyzed.   Finally  for  the  case  of 
Primakoff s  "point  blast"  in  water,  the  solution  is  given 
explicitly  and  graphs  of  various  quantities  are  given  for 
both  the  original  and  the  improveo  solutions. 

It  cannot  be  expected  that  the  range  of  validity  of 
the  "point  blast"  solution  can  be  extended  very  much  by  such 
a  pertTorbation  method.   However  the  qualitative  and  quantita- 
tive nature  of  the  corrections  can  give  an  indication  of  the 
range  at  which  the  solution  begins  to  fail,  and  can  also 
indicate  the  manner  in  which  the  solution  changes  due  to 
atmospheric  pressure. 

Joseph  B.  Keller 
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PERTURBATIONS  ABOUT  STRONG  SPHERICAL  SHOCK  VAVES 

Introduction. 

It  has  been  shown  by  G.  I.  Taylor'  and  others  that  there 
are  solutions  of  the  equations  of  spherical  flow  of  the  form 

u  =  X~'-  ^   U^(r-H) 

p  =  X'2  ^  P  (r-^t) 
t    " 

where  \  is  any  constant  and  U  ,  C^,  P„  satisfy  nonlinear 
■^  o   o    o 

ordinary  differential  equations.   These  solutions,  hoivrever, 
can  represent  a  flox-;  behind  a  shock  only  if  we  can  neglect 
the  pressure  ahead  of  the  shock.   Here  we  shall  study  the 
first  order  effects  of  this  pressure  and  the  modifications 
produced  on  the  original  flow. 

The  problem  of  finding  the  flow  behind  an  expanding 
shock  wave  of  constant  energy  is  reduced  to  solving  some 
ordinary  differential  equations  which  depend  only  on  Yj  '^'^e 
ratio  of  the  specific  heats.   In  the  case  of  y  -  7>  these 
equations  can  be  solved  and  the  perturbations  in  the  flow 
quantities  expressed  explicitly. 

The  difference  betvjeen  the  flov/  behind  a  strong  shock 
and  a  strong  detonation  can  also  be  studied  in  the  same  way. 

Hi®.  Equations  of  Moti on . 

The  equations  for  spherically  symmetric  flow  are: 

u,  +  uu  +  —  p   =0    , 
t     r    p   r       ' 

(rp"Y)^  +  u(pp"^)^  =  0 

G.  I.  Taylor:   The  formation  of  a  blast  wave  bj   a  very 
intense  explosion.   Proceedings  of  the  Royal  Society,  Series  A, 
Vol.  201,  March  19 ?0,  r^.  1^9. 
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where  u  is  the  radial  velocity,  p  is  the  pressure,  p  is  the 
density  and  y  is  the  ratio  of  the  specific  heats  of  the  gas, 
while  r  is  the  distance  from  the  origin  and  t  is  time. 
If  we  introduce  the  variables 

(2)  y^  =  r"^t 

(3)  U  =  Xtr"-^u 

C  =  \tr"-^c  =  \tr"-^  f^ 

•   .    T,  _  -2,2  -2 

P  =  X  t  r   p    , 

where  \  is  any  positive  constant,  we  obtain,  using 

(h)  ^  =   -A  ^  +  ±t 

(4)  ^        ^  ~  +  ~ 

in   (1), 

yi  p 

>2TJ    (1-U)+  UtU^  -r'-^C^  — ^  +  U(aU-l)  +2ax~-^'C^   =   0 

n^M  npv)  *pt        "tc^ 

(5)  -nu^  +__rL(i.Tj)  _2  _''^(i-u)+u-p^-2u  -^-+3aU  =   0 

where    a  =   \"    . 

Equations    (5)    niay  be    rewritten   as 

tP 

Dr^iT^  -  A  -  (i-TT;tu^  -  y'^c^  -^ 

^\  .  ^^t 

where 


1.!  /; 


■         •!     i     / 


a\ 


.^ji-  o 


EiS    ^;^^:'■:J.c•^v'i.M    yo    vj?.;^ 


V  S-.  ^ 


r  ~  ./■• 


;• 


n  w 

..  -g:-". 


•J  ).-• 


-    .'j.-f 


r^'j^ 


A  =  (l-Tj)(l-aTT)Tj  -  {3aTT-2Y~-^(l-a))C^ 


B=  (l-lT)(l-aTT)-(Y-l)[(l-aTj)  -  |(l-a)]TJ-  [  a+y'^l  1-a)  ( l-H)  "^  ]C^ 
(7) 

E=  U-U)[2(1+y)(1-o.tt)  -  3y]  +Y(l-aU)  -  2aC^ 

D=  (1-11)2.  c^ 

There  are  two  Identities  relating  A,  3,  D  and  E, 


(3) 


-A  +  -3Y  B(l-U)  +  D(3aU--3Y  (1-aU))  =  0 

-A  +  Y~"  E(l-IT)  +  D(3aTT--  (l-aUj)  =  0 

Y 


The  special  solutions  for  spherical  flow  known  as 

spherical  "progressing"  v/aves  are  solutions  of  (6)  that  are 

independent  of  t.   They  represent  states  behind  or  ahead  of 

shocks  and  detonations  if  the  shock  or  detonation  front  is 

given  by  )^  =  constant,.  We  represent  any  such  flow  by  U  ( y|^  )  , 

C^ln)  ,  PqCv^)  •   Then  from  (6),  we  find  the  following 

nonlinear  ordlnarv  differential  equations  for  U  ,  C   and  P  , 

^  o'   o      o 

o 


where  A„  =  A(IT^,C  )  etc.. 
o      o'  o 

Prom  the  first  two  equations  of  (6)  v;e  obtain  an 

equation  for  C^  as  a  function  of  U  , 
^  o  o 

dC     3  C 
o  _   o  o 

o  o 

The  solutions  of  this  equation  then  yield,  through  (6), 

F^  and  M  as  functions  of  U  . 
o      *.  o 
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We    shall   consider  here    floirrs   that  are    approximately 
progressing  waves,    that    is   flows  given  by 

(10)  U  =  U^(yj^)    +   eU-j_(Kj_,t)    +   O(e^) 

C   =  Cq(  vx)(l  +  eC-^l  >^,t)  +  0(e    )) 

P   =  p^(  v^)(l+  eP^(  Yl,t)  +  0(e^)) 

where  e  is  some  small  parameter. 

Substituting  (10)  in  (6)  and  equating  first  order  terms 
in  e  we  find 

H^^m  =  (|T  %-^(l|  ^o^i  -  V(i-Uo)t^^it-  y-'c^d^Hp,, 


ncin=(§X^iHl)°  Vi-^^>ii 


(11) 


-1  2 

LV   tC  C^     T 


1-lT     C      2y  l-IJ^   0    It 
o    o      '     o 


n^iA=(^f^  ''i'(§)c  '°'^-"  -^^o't^^it-  (i-^^o)^;'tp,. 


a  A^ 


^--   (Dl=UD^U=U^(yi),  C=C^(VV) 

Conditions  at  a  Shock. 

The  boundary  conditions  on  the  flow  quantities  at  a 
shock  are  given  by 

p(u-z)  =  -  p^z 

2         2 
(12)       p(u-z)  +p  =  p-^z  +p^ 

YP   .  1.    ,2  _  z^      YPi 


i''^)0     '-     (J.   p  ;,-■■?     ->- 


^ 


i  ■; -' 


3^ 


.p  •;     ' 


-  ( :/ ,  _(V      .  '  f    '  n[  p  ) 


.  .4.  '.  ;  ;  o         1 
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x-jhere  z  is  the  velocity  of  the  shock  front,  and  p,  ,  p.,  are 
the  density  and  pressure  ahead  of  the  front. 

In  terms  of  the  variables  U,  c,  P  of  (3)  equations  (12) 
become 

p(U-Z)  =  -  PjZ 

(13)    d(tt-z)^  +  p  =  p^z^  +  \^  yi^\^'^%^ 


2 


(TJ-Z)2  +  (l.^,2^c2  ^  ^222^^^_^2,^2^2^2a^2-2a 


V7here 


(III)     Z  =  Xtr'-'-z 

(15)    [i^   =  (y-i)(y+i)""^ 


2       -1 
^1  "^  '^PlPl 


From  (13)  we  see  that  a  spherical  progressing  wave  can 
represent  the  state  behind  a  shock  for  X  7^  1  only  if  c-.  =  p,  =  0, 
that  is,  only  if  the  shock  is  infinitely  strong.   In  this 
case  the  position  of  the  front  is  given  by 

(17)  vi=H^ 

where  H   is  some  constant  and  thus  from  ilh) 

O  \  r  , 

(18)  Z  =  Xtr"^  ^  =   1 

Vie   assxirae  now  that  the  deviation  in  the  shock  path  is 
small,  and  set 

(19)  Yl=   H^(l  +  eH-j_(t)  )   on  the  shock   . 
Then  from  (lly)    we  find  to  first  order  in  e 

(20)  Z  =  1  -  etli^^ 

Substituting  (20)  in  (I3)  and  setting 

(21)  e  =  ^^c^H^^ 


■■-t- 
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we    obtain 


TT  =    (l-tx^){l  -  etK^^  -er^^"^) 


,^l- 


(22)  C^   =  ^L^(l+|x^)(l-2earH^^+  e      •2'~^^'^_      t^"^^') 

V-    (l+l^    ) 


P   =   p.(l-M.^)(l  -  2etH^^  -  e   -^  t^"^^) 

1+tx 

for     \   =  Hq(1+  eH^(t)).  . 

On  the  other  hand,  expanding  TT(H  ( 1+eH-,  ( t )  ,  t,  e  )  etc.  in 
powers  of  e  we  find 


(23)      U(H^(l+eH^ 


UoO'< 


V^o 


=  ^-^o^'^o 


0(11^(1  +  611-^ 


=  Co(H3 


P(H^(l+eH3^ 


=   ^o(«o 


t)  ),t,e) 


l+eH^(t))j    +   zYi^{Yi^{l+zY.^{t)),t)  +r,{^z'^) 


+  e^'l("^o'^' "^  ^^^--^ 


'^^^oA 


1  o  \^  a  v7  y 


H^(t)  +  0(e'') 


r\-H, 


■ii. 


+  elT-,  (H,,,  t)  +  e 


t)),t,e) 


1  +  eCT(ti    ,t)  +  e 


t)),t,e) 


H-i_(t)  +  C(e'') 


)a=I-I 


D, 


H-^(t)  4-  G(e'') 


n=H^ 


1+  eP-,(K^,t)  +f-f  d; 


l^-o- 


H 


l(t)   +   0(8^) 


■n=H, 


Comparing    (23)    with   (22)    we   have 


(2[|.) 


U^(H^)    -    1    -    M 
C^(H^)    =   1X^(1+  ix-) 


^o^^o^    =   Pl(l-t^    ) 


f  ■ 


■ri 


A    ,  ■'  ''.i^.f   '."j 


.1  .fii.   *-t*i»— 


,:.,-i-    . 


-.11 


liJ  .^-.-i  "   I.,  ; 


,  +  ,f,J',   Arr-^^i 


I.)    "u'j  'Ji 


Pi 


t  - 


i  -■■  I  ^  ^ 


I 


-   "^z 


4.  r  <  '    '  :    -■' 
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Tt  is  then  clear  fron  (9^  that  U  ,  C   and  F  /o-,  are  all 

^  o'   o       O'  '  1 

functions  of  ^/K  ■'.■/hich  defend  only  on  y. 

Frora   the    first    order   terras    in   e ,    xve    find 


Jj^(H^,t)  =  -(p^)  H^(t)-  {l-ii^)tH^j.;t)-  11-/)  t^-^'' 


(25)         C^(H„,t)=-f^)  H^(t)-    tH^^(t)  +  -4^^-t2-2'' 

\    o   /^    „  1+a 

The  perturbation  H^ ( t)  in  the  position  of  the  shock  can 
be  eliminated  and  (25)  reduced  to  tv;o  conditions  on  the  flow 
quantities,  namely  for  V|-  H  , 

"^l^^l  ^  ^1^1  ^  '^1"^    ^   5^t^"^^- 
(26) 

a.2lT^  +  P2C1  +  Y2P1  +  ^21^^-'lt  ■"  ^^2l^^C^t  =  52t2-2':^ 

where  a-,,  (3,,  etc.  are  all  constants. 

To  find  the  flow  behind  a  shock  we  nox-j  have  to  solve 
the  system  of  hyperbolic  differential  equations  (11)  of  third 
order  and  two  conditions  (26)  on  a  space-like  line,  the 
shock  path.   One  more  condition  must  be  prescribed.   For 
example,  we  might  prescribe  the  velocity  of  a  particle 
corresponding  to  a  given  piston  motion  which  does  not  differ 
m.uch  from  that  which  maintains  a  progressing  wave.    This 
would  involve  extensive  com'outations  although  it  can  be 
reduced  to  a  second  :">rder  problem. 

Alternately  and  more  naturally,  wo  may  prescribe  the 
total  energy  contained  in  the  shock  wave,  that  is,  the 
energy  imparted  at  t  =  C.   This  is  essentially  equivalent 
to  prescribing  one  condition  on  the  t-axis,  namely  tho.t  no 
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energy  is  added.   From  the  general  theory  of  hyperbolic 
eqiaations  we  will  have  the  correct  number  of  conditions  to 
determine  the  i^roblcm.   In  fact  this  problem  can  be  solved 
in  terms  of  the  solutions  of  ordinary  differential  equations. 
In  the  case  of  a  shock  wave  in  water  it  can  be  solved 
explicitly . 

Shock  Wave  of  Const  an  t  ^  En  c_r  g^  • 

First  wo  note  that  equations  (11)  with  the  boundary 
conditions  (25)  have  the  special  solutions, 

%  =  t2-2^()/^^(  v^)  +  hXi2(n)) 


(27) 


C-L  =  t2"^''-(X2lK)  +  hX22(  VV)) 

H^  =  ht'--2a 

where  h  is  an  arbitrary  constant  and  X.n  i  '  )Cr>-\)     Xqi  -'-^  ^^^ 
solution  of 


2   1-n,  ;^^.a^££!^,l.,;^. 


l-"o   =0 


dX 


n-dT^  =(l)°Xi  +  (§)°  C„X2  -  2Y{l-a)D-lXl 

which  satisfies  the  initial  conditions 


I  if 


A 


'  (• 


'  I  ■'  \ 


K  -, 


!  '■  '  ^ 


9. 

(29)  Xai^^^o^  =  -¥^^ 

1+11 

and  Xip'     ^?P'     "^TP  -^^  ^''^^    solution  of  (28)  which  satisfies 
the  initial  conditions 


>^12(^o'  =  ""(d;')     "  2(l-tx^)(l-a) 


■A 
o  1 

r,  ; 


o 


(30) 


>^22(H^)  =  -(?)     -  2(l-a) 


n=H„ 


'■^o 


^32(^^0^  =-   D       -  ^+(1-^^   • 


A  method  for  reducing  the  third  order  systom  (26)  to 
one  of  second  order  is  contained  in  the  Appendix. 

It  turns  out  that  the  so?.ution  (27)  is  just  the  solution 
which  satisfies  the  condition  of  constant  energy  if  h  is 
chosen  appropriately. 

Let  the  total  energy  contained  in  a  spherical  shock 
wave  at  any  time  be  E(t)  where 

(31)        E(t)  =  Utt  r    U  pu^  +  -— y)  ^"^  dr 

o 

Here  R(t)  is  the  position  of  the  shock.   In  terms  of  the 
variables  (2)  and  (3),  (3I)  becomes 

H(t) 

00 


(32)    E,t,  =-W  J  ''■  (I  in^  .  ^)  P  ^-Sa-1  .5-2  „ 


3  Y- 
where   y^=  H(t)  represents  the  shock. 


h^IlE^it) 
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In  the  case  of  an  infinitely  strong  shock,  uaevc   x)     =   0, 
E(t)  vjill  be  constant  only  if  a  =  2/5.   rnen 


(33) 


E  =  E. 


Mt)  2 

CO        o 


12^ 


P^v|"3  ^Yl 


If  we  now  consider  that  second  order  terms  in  p.,  can  be 
neglected  and  use  the  approximations  (10)  and  (19)  wc  obtain 
H, 
(3^-)   E  =  -  g^ 


n2 


?F   J     \?  ~?  "^  y-T]  ^o  ^      "^^^ 


~     -i  -)"r'     KW  ^W  -I 

12S         O    1 


Po  n 


-3 


J^-'^o 


32Tre 
12^ 


o 

ilTtp. 


o         3  (  Y"  1 .'       o 


From  the  first  order  terms  in  e  using,  from  (21), 


=  ^2   2  .^2a  ^  25 


■1  u^/3' 


e  =  X   c-L  H^-  =  '^-   YPiPi^  H- 


we  obtain 


(35)  k;\ 


'i  yUo 

2  — ? 


y-l  /   o 


Jn-H, 


H 


00 


^   -'o       ^ 


■^^—^■•Wi^l'  ^"^  dV^ 


2  YH'-l)        o 


or,  using  (2I4.)  , 
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(36)    Hj  -'■^r^  + 


i-:^*^)'! 


00 


=0   "" 


yu 


.-3 


|!V'wKi(^p(t)--l7(7^-'/^ 


We  can  now  show  that  (27)  is  just  the  solution  of  (11) 
and  (25)  v/hich  satisfies  (3  6)  if  h  is  chosen  appro-oriately . 

Substituting  (27)  in  (36)  using  r~  t  =  H^  to  first  order 
on  the  shock  yields 


ht 


6/5 


2(1-M-^) 


Y 


:T" 


00 1^0       ^ 


+ 


■■jfa.-^x...Hiftr"'ft)  -"^'I  ( 


I'-r    0-2 
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Xu 


N-21/5 


o 
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(37) 
vjhc  re 


A3_h   +   A2    =  A3 


00  L  o 


Yll 
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x..-^M  'fMT  <^) 
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(38)      A. 


iS^-ilXoi-^Xn-^oX 


CD 


y 


1    ■••-^o  / 


v-21/5 


JL 


H 


o 


A      =   - 


7   Y(Y-l) 


1  --■■  ) 


.u  \ 


.  u      Of,  ;    ■>!     1 


J  •-•■,' 


/    •• 
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J'  s 
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12. 

Equation  (37)  can  always  be  solved  for  h  provided  A-,  i"   0. 
Thus  the  first  order  iacrerient  in  the  energy  vanishes  and 
we  have  a  complete  solution  of  the  problem  in  terms  of  "^  .. 
etc . 

From  (7),  (9)  and  (2ii)  we  see  that  IT   and  C   and  P  /p, 
are  functions  of  y^/  H   only  and  independent  of  p,  while 
F/p,  is  a  function  of  V|/h  . 

Thus  the  coefficients  of  the  equations  (28)  are  functions 
of  V^/H^  and  hence  by  (28),  (29)  and  (30)  the  functions  Xll ' 

X21'  X31'  Xl2'  X22'  X32  ^^P-"^  °^  ^A'o  ^''ly- 

Thus  A^  and  k^   are  constants  depending  only  on  the 

solutions  of  fixed  differential  equations  with  fixed  initial 
conditions,  depending  only  on  y 

Substituting  (2?)  in  (10)  and  tnen  in  (3)  we  find 


u 


2  £  ,„  ..t6/S, 


=  =  I  t-3/5  vi-^/5,^,|  ,,3A  >i:2/5c^(X,,.hX22) 


Note  that  for  fixed  Y\^  the  perturbations  in  u  and  c 
increase  with  r  vjhile  the  pcrtiirbation  in  p  is  constant. 

Novj  the  pressure  dies  out  behind  the  shock  and  the 
maximum  pressure  occurs  at  the  shock;  thus  from  (22)  and  (3) 
we  find 

2 

Tn  other  words,  the  maxim.ijm  pressure  is  the  pressure  of  the 
Taylor  point  blast  wave  plus  a  constant. 
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13. 

Sh o_c kxja V e    of   Constant   Ene rg v   In  Wa t e r  . 

In  the  case  of  an  explosion  in  water  the  functions  X-i-, 
etc,  can  be  calculated  explicitly. 

Here  we  have  F  =  A{ (-^)   -1)  where  A  =  30C0  atmospheres 

and  Y  =  7'   Then  we  replace  the  last  equation  of  (3)  by 

2  2-2 

P  =  \  t  r   (p+A).   In  this  case, 

-   .  -2  /^ 
(39)   F^=  1-lx-  ,  cl=   /{l^/)     ,   P^=  p^n-li^)    (^^) 

is  the  solution  of  (9)  and  (2[|.)".   Then  A   =  0,  B   =0  v;hile 
E  and  D   are  constants. 

The  differential  equations  (28)  then  have  constant 
coefficients  in  the  lovjest  order  terms  and  the  solutions  are 
powers  of  r|_  .   Satisfying  equations  (29)  and  (30)  we  find, 

:.i;738  .       ^0.3212 


/'        > '5, 4- coo  ^_    xL 

Xii  =  -1.00197  i^^J  -   .5605  (^) 


•1.6000  ^      N5.i!-788 


X21  =  2.2689  (f-1  -5.4262  l^] 

v"  o  ^  V    o  y 


.'.3212 


.  2.3597  (^) 

/  «   \-l.6'"00  /  „  ,5.14.728 

X,i  =  -1.5126     f  -2.3^71  (f) 

■v      0/  V    o  ' 

/y,    f.3212 
+  1.1311  (  ^  ) 
V'^o/ 

/vi  n5. 1+738  /^  ^0.3212 

X,2  =  -0.3885  (1^)  +  0.0885  (jL) 

/^,  ^-1.6^oo  /.^  •v5..'+738  /- ^   .0.3212 

X22  =  1-3553  (I-)  -2.1039  [^)  -  -45ll|(^) 

•1.6- vO  ,  ^  .5.4783 


X23  — 0.9035  f^^)"    •   •■•    -0.9100  (^j.^j' 


-   0.1865  (^) 


0.3212 


This  special  solution  Is  due  to  H.  Primakoff , 
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Substituting  (39)  and  ([lO)  in  (38)  and  tbon  (37)  we  find 
t^iat 

h  =  -3J.1.938 

and    finally,    from   (27), 
(^1)      n^  =  t°/^i.3553(7f)  -.3697  (jf-)  > 

■\{  /w  N-^-.600(  /      ^.1^738  /^   xO.32121 

C-L  =   r  J  -1.61^2    (JLj  +.8323(J-j  +1.8326^^]  j 

Pi   =  r^J  -2.1,662  (^|Lj  +i.92[,2(|Ly  +L|..l:.369(f  ■)  V 

H^  -  -3.i]-938  t^/-''        . 

In  Figures  1,  2  and.  3  wo  have  plotted  the  path  of  the 
shock  as  a  function  of  tirae  and  the  maximum  pressure  as  a 
function  of  time  and  distance. 


^ '^ ^.  Behind  a  Do tona t i on . 

The  third  shock  condition  in  (12)  miist  bo  modified  to 
include  the  chomical  energy  of  the  cctoaation.   VJo  then  have 
for  the  conditions  across  the  front 

(1|.2)       p(u-z)  =  -  p-|^z 

2         2 

p(u-z)  +  p  =  p  z  -;-  p^ 

n       p   —     2      YP-i       — 
P(y-1)    2'  2    p^(y^-1;     1 

vjhcre  E  and  E,  are  the  energy  of  formation  per  unit  mass  of 
the  burnt  and  unburnt  material  respectively,  and  y..  is  the 
ratio  of  the  specific  heats  in  the  unburnt  gas.   In  terms  of 
the  variables  (3)  and  (ll^.)  the  last  condition  of  (I4.2)  becomes 
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ti^(TJ-Z)^  +  (1-1J,^)C^  -  n^z^ + 


15. 

2 


[J- 


o 


^2   2x-2^2 


If  vie   DertiTrb  about  a  stronp:  detonation,  i.e.  setting 
2 

^(  1-jj,  )C  +  E  -  E-,  =  0  we  can  find  the  undisturbed  flov;  as  a 

^^o 

spherical  wave  and  the  perturbed  flow  satisfies  (25)  with 

2x-2 
different  constants  for  the  coefficients  of  r 

In  this  case  there  are  again  special  solutions  of  the 
form  (27)  where  now  the  initial  conditions  (29)  must  be 
adjusted  appropriately,   Hov;ever,  in  this  case  we  need  one 
more  condition  on  the  flow  and  the  special  solution  will  not 
in  general  satisfy  it. 

For  Y  =  2  we  have  U  (H  )  =  "^o^^o^  or  u  =  c  at  the  shock. 

The  unperturbed  detonation  is  then  a.  Chapman- J ouguet 

detonation.   If  the  pertitrbed  tlo\-j   is  also  behind  a  Chapman- 

Jouguet  detonation  we  obtain  a  relation  between  Ut  (H  ),  0^(1-1  ) 

i^  o  '   1^  o 

and  H, (r).   This  relation  can  be  satisfied  by  solutions  of 
the  form  (27)  for  an  appropriate  choice  of  h. 


Appendix. 

In  this  Appendix  we  will  show  hov:  equations  (26)  can  be 
reduced  to  one  second  order  equation.   In  general,  equations 
(6)  and  the  shock  conditions  (25)  can  be  reduced  to  a  second 
order  equation  Involving  the  unknovm  shock  function  H-,(x). 

Prom  equations  (o)  we  see  that 
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and  from  (11)  then. 
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1^2   /B^°    _„l_rAN°W     _iL.  (1-c^)  c"^V 


-1  ^x^ 


l-U 


n 


d 


=  [r-  (31  -  ^  (^D  ^- 


•1  2(l-a)   ^ 


1-TT. 


«>^i 
'1^^- 


((5)°^2(l-u„)(l.a)E;l)Xi 


If  WO  introduce  new  dependent  variables. 


(A. 3)        ?i(  Vt)  =-  Xi(  A)A^o(l-^^o^ 


2 


^2^^!^'  =7-T  X2^^^  ^  ^i(vv)u, 


-1 


C3(va)  =  y"   X3^^)  +  ^i(VV)TJ, 


equations  (A. 2)  reduce  to. 


(^'^)   (l-^^o^^dl^ 


=  W. 


2C' 


o 


2U^-1  A    1-IJ   ^,  xo     ,  C   /.>o     / 
o  o  ^ 


\0  C 


2C 


/^AY  ^o  Y-1  ^   ,   """^o  ,,   .^ 


o  o 


.1  -- 
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^-^^o  .     ^'^2   _    ,    .    .      2 


1-U    d^  . 

ir^  1  dlf  "  '"^  "*■  2y"  (A-l)  +  2(l-a))^^  -u;-^(2-2g)^3 

By  introducin.'^  new  dependent  variables  we  can  reduce  (A.Ii) 
to  a  much  simpler  form,   V7c  set 


where  Hw.  is  any  fixed  value  of  r^.      Then  (A.I4.)  reduces  to 

(A. 6)      1^0  4y-^  Ko^i-^V.2^2(l-a)  ^  C3  =  0 

o 

where  K  and  L  are  functions  of  y, 

{2U^-1     1-U     ^,  xo 


-^-5i^=o(§|}-^=-)^ 


o 


and  ^ip  and  \x^   are  constants, 


ii-f 


1.4-1.  t'. 
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(A. 8)        1x2  =  3  -  2^'-^^  ^   2y-1  -  2{l-a) 

^^3  =  3  -:^A.:^-  2(l-a)    . 

The  solutions  of  tho  third  order  system  (A. 6)  can  also 
bo  expressed  in  terms  of  the  solutions  of  a  second  order 
differential  equation.   Wc  introduce  the  function  G  x^here 

(A. 9)  ^  ^  ^  h 

(A. 10)  G(Y^)  =  0 

where  Y  is  some  fixed  value  of  y. 

Then  from  the  last  equations  (A. 6)  wc  obtain 

(A. 11)  ^2^y'  "  ^2^Yo^  "-   ^2^^^^ 

(A. 12)  K^ij)  =  ^3(^0)  -^  t^3c-(y)     • 

Substituting  (A. 11)  and  (A. 12)  in  the  first  equation  of  (A. 6) 

yields 

2 

(A.13)   D^  ^  ^  K^  §  +  gU^L^   +  (2-2a).^3  ;^  ) 
dy        0^     \  Uq 

o 

NOW  £,^   and  f.^  arc  prescribed  for  y  =  Y^  by  (29)  or  (30) 
Thus  wo  have  a  differential  equation  (A.13)  for  G  and  two 
bound  iry  conditions  (A.9)  -nd  (A. 10)  at  y  =  Y^. 
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